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Abstract 

We present the results of a new investigation aimed to identify the pulsational 
characteristics of the Southern (3 Cephei star ALS 3721. Spectroscopic and 
multicolour photometric data were acquired at the South African Astronom- 
ical Observatory (SAAO), South Africa in 2011. Frequency analysis showed 
that the oscillations of ALS 3721 could be attributed to the two main fre- 
quencies with higher significance. Stellar parameters and projected rotational 
velocity obtained by the spectra were also used to perform photometric mode 
identification. In order to determine spherical harmonic degrees (7), a princi- 
pal method was followed by comparing the observed light amplitude ratios in 
different passbands with those computed from non-adiabatic pulsation mod- 
els. In general case, therefore, the spherical harmonic degrees corresponding 
of the frequencies were found in the expected /3 Cephei range. 
Keywords: stars: oscillations (including pulsations) - stars: early-type - 
stars: individual (ALS 3721) 
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1. Introduction 



The (3 Cephei stars are defined as a group of early B type near main 
sequence pulsators (B0.5-B 2.5) with the masses between about 8 and 18 M Q 



( jStankov fc Handler! . 120051 ). They show both mono and multiperiodic short 



term pulsational variability in the periods of several hours. Although their 
general pulsational behaviour is considered by low order p (pressure) modes, 
a number of them pulsate in low order g (gravity) modes and pulsations are 
driven by k mechanism throug h iron group elements in the partial ionization 



layer s of their stellar structure (ICox et al. 



1993 



1992 



Dziembowski. Moskalik. fc Pamyatnykh 



Dziembowski fc Pamyatnykhl . 



1993T ) . Like most other type 



of pulsating stars, /3 Cephei stars are represented with a well known instabil- 
ity strip in the Hertzprung-Russel Diagram (HRD). However, the blue part 
of this instability domain is not well populated since there is lack of chance to 
detect unstable modes for such massive (> 25 M ) stars with strong stellar 
winds (jPamyatnykhl . Il999[ ) . In a similar way, it is also very difficult to observe 
very low amplitude modes at the red part of the strip from the ground based 
observations. A detailed revie w of the Kepl e r obse rvations of B type vari- 



ables were recently reported by 



Balona et al. 



(120 111 ). They pointed out that 



instability domain among B type pulsating stars may also be divided into 
two subgroups in the HRD as the metal opacity bump occurred in different 
layers for each star group: i) high-luminosity Cephei and ii) low-luminosity 
SPB (Slowly Pulsating B stars). In this case, the most recent asteroseismic 
studies for these type of variables are being more important to understand 
their strange physical structures. 

The star ALS 3721 is a member of the Southern open cluster NGC 6200 lo- 
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cated in the constellation Ara. The first observations of the cluster was made 
by 



Fitzgerald. Jackson, fc Moffat 



( 119771 ). The authors derived the colour ex- 



cess E(B-V)= m 63 ± m 07 and the distance modulus ll m .9 ± m .2 of the 
cluste r. ALS 3721, o n the other hand, was classified as B0 .5 

(119771 ). 



((19631) and B1III by Fitzgerald. Jackson, fc Moffat 



II by 


Whiteoak 


Pigulski 


( 


2005) 



indicated that the s tar has Cephei type oscillat ion characteristics based on 



the ASAS-3 survey. 



Pigulski &: Pojmahskil (120081 ) analyzed the ASAS-3 data 



in order to determine the oscillation frequencies of ALS 3721, and they re- 
ported the three independent fr e quenc ies with highest amplitude. Recently, 



Ulusoy. Engelbrecht. fc Giilmezl ( 120121 ) analyzed the new photometric data 



taken from the SAAO by using the Lomb-Scargle periodograms and they 
referred to ten oscillation frequencies that are reached in the B light curve. 

In this paper, we present a detailed study of the pulsation properties 
of ALS 3721. Spectroscopic observations and analysis of the target are ex- 
plained in the next section. In the third section the details of the photometric 
observations, reduction and frequency analysis are given while the mode iden- 
tification for each detected frequency, which is the main goal of this study, 
and its results are stated in the fourth section. The study is summarized and 
results are discussed in the last section. 



2. Spectroscopy 

The spectroscopic data were used to derive the fundamental stellar at- 
mospheric parameters and the rotational velocity, vsmi. The observations 
were carried out on 18 June 2011 with the Robert Stobbie Spectrograph 
(RSS) longslit instrument mounted at the Southern African Large Telescope 
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(SALT), Sutherland, South Africa. The specifications for the RSS instru- 
ment yielded R ~ 3000 spectral resolution of given characteristics related to 
the spectra taken for ALS 3721. Four spectra were obtained with exposure 
times set to 60 s and signal-to-noise ratios were achieved from 80 to 100 with 
a spectral range of 3925 and 5990 A. Reductions for the RSS spectra were 
performed using the IRAF-0NEDSPEC and the TW0DSPEC packages, also the 
normalization of spectra were carried out by means of the C0NTINIUUM pack- 
age. In order to obtain atmospheric parameters of ALS 3721 we compared 
the low-resolution averaged SALT spectrum with theoret ical spectra from 
the TLUSTY/SYNSPEC codes given in the BSTAR2006 grid (ILanz fc Hubenvl . 



20071 ). The BSTAR2006 grid consists of 1540 metal line-blanketed, non-LTE, 



plane-parallel, hydrostatic model atmospheres. The synthetic spectra were 
calculated for effective temperatures between 15000 and 30000 K with a step 
equal to 1000 K and surface gravities ranging from 1.75 to 4.75 dex with a 
step equal to 0.25 dex. In our fitting we adopted fluxes calculated for solar 
metallicity. The spectra are available for microturbulence velocity equal to 2 
kms -1 . The model spectra were convolved with rotational and instrumental 
profiles. 

The determination of surface gravity log g and effective temperature T e g 
was based on the whole available spectrum. The important diagnostic fea- 
tures are clearly visible. Four Balmer lines (He, H<5, H7 and H/3), He I lines 
(4026, 4388, 4471, 4713, 4921 and 5876 A), silicon lines Si III (4552, 4568, 
4575 A) and blended lines of other elements, like Si II, O II, C II, N II, S II and 
FeUI were identified. The best fit of the model fluxes to the observed lines 
was obtained as a result of minimizing residual standard deviation of the fit 
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Figure 1: Observed (black) and synthetic (red-fitted) spectra of ALS 372 f for selected 
region. The observed lines fitted to model fluxes with T D g — 24000 ± 1000 K and log g = 
3.5 ±0.1. 

(Fig. [[J. The analysis resulted in T eS = 24000 ± 1000 K and logg = 3.5±0.1. 
Additionally, the projected rotational velocity v smi, was derived to be equal 
to 100 ± 20kms _1 and it is obvious that higher microturbulence velocity is 
needed. 



3. Photometry 



The multicolour photometric data were used to calculate the amplitude- 



colour relation for the two stronge st frequencies of ALS 3721. 



achieved using the data taken by lUlusoy. Engelbrecht. fc Giilmea (120121 ). 



his was 



They carried out the observations with the 75-cm and the 100-cm telescopes 
at the Sutherland station of the SAAO. Since the data obtained from the 
75-cm telescope is scattered than expected and a large time gap between two 
observing runs from the different instruments exist we included only data 
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from the 100-cm Elizabeth Cassegr ain telescope to our investigation. Th e 



details of the study are explained by 



Ulusov. Engelbrecht. & Gulmezl fl2012[ ) 



During the observations, a total of 780, 866, 871, and 853 observing points in 
the Johnson U, B, V and / filters were collected, respectively. Data reduction 
was performed following standard IRAF routines. In the pre-reduction pro- 
cess, the bias and the dark frames were subtracted, the flat-field correction 
was applied to each of CCD images. Instrumental magnitudes for all stars 



in the field of view were measured by aperture photometry techniqu e. A 



of the images were processed using the DA0PH0T II package (iStetsonl . 119871 ) 
which is also included in IRAF following the standard tasks of the package. 
The data were reduced relative to only one comparison star, GSC 8330 2327, 
due to showing the least variability in the observed CCD field. The observed 
light curves of the star in the UBVI filters can be seen in Fig. [2J 



Oscillation frequencies were computed with the PERI0D04 (ILenz fc Breger 



20051 ) software for each filter. The sofware package is based on the classi- 
cal Fourier analysis and applies least-squares tecnique to input data by the 
fitting formula: 



(1) 



where Z is the zero level of magnitudes, A4 is the amplitude of corresponding 
frequency, Qi is the frequency value, t is time and $j is the phase of given 
frequency. The si gnal-to- n oise r atio S/N > 4 threshold was adopted as a 



criterion following 



Bregerl ( 120001 ) . The theoretical frequency resolution was 



found to be 1/T = 0.4494. Since (3 Ceph ei type oscilla t ion fr equencies are 



generally detected in the range of 4-20 d 1 (IBalona et al.l . 1201 ll ) we restricted 
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our analysis on the frequency range between d^ 1 to 15 d _1 . 

The analysis completed with the detection of the independent and com- 
bination frequencies. The results correspond that two frequencies with the 
highest amplitudes can be assigned for the oscillation of ALS 3721, /i=4.9457, 
/ 2 =5.3928 with the amplitudes of A 1= 0.08404, A 2 =0.02517, respectively 
The amplitude spectra of 5-filter data before prewhitening of any frequency 
between mentioned interval is shown in the middle panel of Fig. |3j Peaks 
with higher amplitude around 5 d _1 is clearly seen from the figure. The 
bottom panel represents the final spectra after overall process. The obtained 
frequencies, amplitudes, phases and S/N values for all light curves after the 
analysis are sorted in Table [TJ Signal-to-noise ratios were computed in the 
interval of 2 d" 1 . Light variation of the star in the UBVI filters correspond- 
ing of multiperiodic oscillations is shown for six different nights in Fig. |2j 
Solid lines in the figure show the agreement between analysis results and the 
observational data. 

We conclude that ALS 3721 oscillates in t he frequency range of Cep 



type oscillators which were also confirmed by 



Pigulski fc Poimahskil feOOgJ ) 



previously using ASAS-3 data. Tab 
light curve to that of determined by 



e [2] compares our ana l ysis re sults of V 



Pigulski fc Poimahskil ( 120081 ) . 



4. Mode Identification 

There are different mode identification methods that can be used to iden- 
tify the pulsation mode geometry of a star. An oscillation mode is specified 
by three spherical harmonic quantum numbers; I the degree of the pulsation 
mode, m the azimuthal order, and n the radial overtone. In the photometric 
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Table 1: Computed frequencies, amplitudes and phase shifts during the solution. Signal- 
to-noise ratio per frequency is given in the last column. The standard deviations in the 



last digits are given in parentheses. 

ID / (d' 1 ) A (mag) S/N 

u 



Si 




4.94568(25) 


0.12067(124) 


0.8872(16) 


65. 


Si 




0.0719(13) 


0.02408(124) 


0.7860(82) 


10. 


h 




5.3951(10) 


0.02965(124) 


0.9077(67) 


14. 


Si » 2/! 




9.8727(12) 


0.02616(124) 


0.2369(77) 


14. 


h 




10.4076(24) 


0.01216(124) 


0.9613(157) 


7. 


fa 
J D 




3.9682(10) 


0.01118(124) 


0.9042(166) 


5. 


St 




11.0997(33) 


0.00851(124) 
B 


0.3956(216) 


4. 


Si 




4.94565(8) 


0.08404(55) 


0.3131(5) 


54. 


Si 




5.3928(6) 


0.02517(55) 


0.1529(36) 


15. 


S3 




0.0719(9) 


0.01521(55) 


0.9118(60) 


20. 


Si W S3 + 2/2 




10.8750(6) 


0.02039(55) 


0.9725(44) 


24. 


S5 ~ Si + Si - 


/2 


10.4143(11) 


0.01206(55) 


0.0711(33) 


12. 


Se « Si 




4.9726(1) 


0.01236(55) 


0.4157(9) 


8. 


Sr 




13.8882(13) 


0.01095(55) 


0.9684(82) 


4. 


Ss 




1.3707(15) 


0.00806(55) 


0.4107(99) 


8. 


So ~ a/2 - Si 




11.2030(10) 


0.00826(55) 
V 


0.3873(102) 


10. 


Si 




4.94563(14) 


0.07809(44) 


0.0632(9) 


85. 


S2 




0.02472(11) 


0.02607(44) 


0.1498(7) 


18. 


S3 




5.3905(5) 


0.02224(44) 


0.7944(32) 


23. 


Si » 2/! 




9.8706(6) 


0.01927(44) 


0.5048(37) 


21. 


S5 w Si + S3 + 


2/2 


10.4121(12) 


0.00929(44) 


0.3159(77) 


11. 


Se w /1 




4.97259(9) 


0.01297(44) 


0.3456(6) 


14. 


/7 « 3/ 6 




14.8993(10) 


0.01055(44) 


0.2958(69) 


8. 


/s 




1.3167(15) 


0.00657(44) 


0.9132(94) 


5. 


/9 




11.1514(14) 


0.00791(44) 


0.3691(90) 


10. 


S10 ~ Si + 2/4 


"/5 


14.2365(19) 


0.00531(44) 


0.3175(121) 


4. 


Su 




3.5413(34) 


0.00411(44) 
I 


0.9863(220) 


4. 


Si 




4.94338(22) 


0.05811(53) 


0.8090(14) 


51. 


Si 




5.39502(54) 


0.02422(53) 


0.7542(35) 


19. 


S3 




0.0719(10) 


0.01366(53) 


0.7355(62) 


5. 


Si » /3 + 2/2 




10.8727(7) 


0.01924(53) 


0.1160(44) 


9. 


/*5 « /l + /4 — 


/2 


10.4143(10) 


0.01209(53) 


0.2394(70) 


6. 


/6 




4.1592(13) 


0.01001(53) 


0.3548(84) 


6. 
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Figure 2: Observed (dotted) and theoretical (solid) UBVI light curves of ALS 3721 for 

9 

six different observing nights. 




-15 -10 -5 5 10 15 




Frequency (c/d) 

Figure 3: The frequency spectra of ALS 3721 in the B filter. The top panel shows the 
spectral window. The middle panel represents the periodogram yielded before the analysis 
while the bottom one shows after pre-whitened all listed frequencies in Table [TJ 
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Table 2: Comparison of Four ier analysis results of V light curve to PP08 which stands for 



Pigulski fc Poimahskil (|2008l). 



ID This study PP08 

/i(c/d) 4.9456 4.9488 

Ai(mag) 0.0781 0.0872 

/ 2 (c/d) 4.9246 
A 2 (mag) 0.0151 

/ 3 (c/d) 5.3905 5.3989 

^(mag) 0.0261 0.0116 

ft ~ h + /a(c/d) 9.8706 9.8734 

A 3 (mag) 0.0193 0.0105 

* The combination frequency contains / 2 which we could 

not obtain but published by PP08 previously. 

case, spherical harmonic degree, I can be derived by comparing observed am- 
plitude ratios, or phase differences for a given passband with the predicted 
ones. Hence, mode identification from multicolour photometry is an essen- 
tial requirement for modelling the oscillation s in pulsati ng stars. For our 



study, we used the FAMIAS software package ( jZimal 120081 ) to determine the 
wavelength dependence of the photometric amplitudes corresponding to the 
mode degrees of the first two strongest frequencies found in ALS 3721. In 
FAMIAS, the photometric module applies the aforementioned method with 
precomputed model grids of non-adiabatic stellar models for a certain range 
of stellar parameters in a number of photometric passbands including the 
Geneva, Johnson/Cousins and Stromgren systems. Whereas the determina- 
tion of the I degrees strongly depends on pulsational input parameters, our 
calculations are performed in a wide range of the model grids to fit the best 
estimation of the I degrees. 
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For ALS 3721, we calculate the theoretical amplitude ratios of the < I < 
3 modes for the models assuming a mass range between 10M Q < M < 14M Q 
and effective temperatures in the range of 23000 < T c g < 25000 K. The 
derivation of the theoretical values of different I degrees in a certain range 
of T eS and log a were obtained by using 1.6 < M Q < 20, ZAMS to TAMS 
and A04, OPAlH The theoretical expression of the model (with the assump- 
tion of zero-rotation and static plane-p arallel atmospheres) included in the 



software package was represented by (IDaszvriska-Daszkiewicz et al. 



Daszyhska-Daszkiewicz. Dziembowski.&: Pamyatnykh 



2002 ; 



20031 ) . The precom- 



puted grids of the atmospheric parameters have been calculated by applying 
Kuruc z and NEMO a tmospheres and the quality of the solution is represented 



by x 2 dZima 



2008). As shown in Fig. HJ the amplitude ratios were plotted as 



a function of the central wavelength normalized to the U filter for the highest 
amplitude frequencies /i=4.9457 and /2=5.3928. For the primary pulsation 
frequency, fi, the relation between the predicted values and the observed 
ones are not too much distinctive for the masses of 10M Q and 14M . It 
may be concluded that the amplitude ratios are intermediate between I = 
and 1 = 1. In contrary, the masses between the values of 1OM and 14M Q are 
remarkably represented with the curve of / = 0. The calculated amplitude 
ratios for f2 except for the mass of 14M are matched with the observations 
and lead to value of / = 1 but the results are not so clear as those of masses 
for 14M . Consequently, mode identification is inconclusive for this mass 
and it could be considered between I = and 1 = 1. 



1 0PAL opacities computed with non-adiabatic Warsaw- New Jersey/Dziembowski code 



by J. Daszynska-Daszkiewicz and P. Walczak (http://helas.astro.uni.wroc.pl). 
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Figure 4: Observed and theoretical UBVI amplitude ratios for ALS 3721. The colour 
amplitudes are given as a function of the central wavelength normalized at the the U filter 
for fi (left) and ji (right). The coloured lines in the panels represent different spherical 
harmonic degrees corresponding to different stellar models in the range of 23000 < T e g < 
25000 K and T off , and logg 3.2 < logg < 3.8. The red (solid) lines indicate 1 = 0, the 
green (long-dashed) 1 = 1, the dark blue (dashed) I = 2 the purple (dotted) I = 3 . The 
models computed a mass range of between 9 and 15 Mq. The black crosses represent the 
observed values with their lsigma standard^ deviations for fx and f%- 



5. Discussion 

ALS 3721 was confirmed to be possessed (3 Cephei type pulsation char- 
acteristics. Frequency analysis based on photometric observations resulted 
with the two reasona 
previously proved by Piuulski <k' Pojnuuiski (2008). Spectroscopic observa- 



jle frequencies, /t =4. 9 457 a nd / 2 =5.3928, which are 



tions and analysis indicate that ALS 3721 has a high rotational velocity ( 
vsini=100 ± 20kms~ 1 ) with derived effective temperature T e g- = 24000 ± 
1000 K and log g = 3.5 ± 0.1. Additionally, the luminosity of the star is 
found to be log7^=3.22 based on the calculation from the distance modulus 
of NGC 6200, m- M = 13 m 36, given by WEBD^H database. During the cal- 



culation, the values of appare nt magnitude of the star and the s olar absolut e 



magnitude was adopted from 



Pigulski fc Pojmahskil (120081 ) and 



Coxl (I2000h . 



respectively. ALS 3721 is a member of the Southern open cluster NGC 6200 
and it is an early age cluster. In this case, if all results were taken into con- 
sideration it might be concluded that the star was located on the ZAM S and 
still at the stage of core Hydrogen burning in the instability strip. 

One of the bottlenecks for existing mode identification methods is to 
determine the order of spherical quantum numbers in an accurate way. Re - 



cently, with new generation space missio ns like Kepler { 



CoRo T (IBaglin et all . |2006J) and MOST (IWalker et al. 



Borucki et al. 



20101), 



20031 ) stellar variabil- 



ity is re-identified by analyzing unprecedented data sets which makes possi- 
ble the detection of the very low amplitude oscillations with a high precision 
level. Alas, current space based data are not compatible with existing models 



1 http:/ / www. univie. ac. at/webda/ webda.html 
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including both spectroscopic and photometric methods and thus astreoseis- 
mic studies are still in need of the ground based data. From the ground 
based observations, /3 Cephei stars are distinguished with typical frequencies 
in the high frequency domain. Considering that, these stars are relatively 
more convenient targets in order to test the reliability of the stellar models. 
Since it is assumed that they have no deep surface convection zones, several 
physical conditions would be negligible in the model computations. Because 
of this, seismic interpretation of the pulsation modes can be more possible 
for at least a few frequencies, however, this is not a satisfying solution in 
some cases such as rapid rotation. What is more, model computations are 
more complicated than expected in the presence of rapid rotation. In other 
words, if the star rotates very fast, the effects of rotation on eigenfrequencies 
should be additionally considered for a trustworthy seismic modelling. For 
the j3 Cephei stars , another importan t effect on the detection limits of th e 



oscillation modes were also reported by 



Daszyhska-Daszkiewicz et al. 



J2002I). 



They have found that the determination of higher I modes (I > 2) only from 
multicolour photometric data may not be expected. Consequently, such I 
modes may be expectable in highly densed frequency spectra obtained with 
space based instrumentation. On the other hand, the detection of intermedi- 
ate I modes (/ < 2) are more likely possible by using spectroscopic line profile 
variations and photometric amplitude ratios. In our study, the spherical har- 
monic degrees have been computed for the two of strongest frequencies of 
ALS 3721. In addition to the photometric data, spectroscopic data were also 
achieved to extract fundamental stellar param eters. For mode identification, 



we used FAMIAS software package (jZimal . 



20081 ) to obtain amplitude ratios cor- 
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responding to I degrees assuming the masses between 1OM < M < 14M . 
Model outcome for the mass of 1OM and 14M as seen in Fig. HI the dis- 
tinction between I = and I = 1 modes for f\ is not clear in the amplitude 
ratio. The results for f 2 the amplitude ratios are consistent for a mass range 
of between 1OM < M < 13M and reproduced by the curve for 1 = 1. 
Nevertheless, observed and calculated colour amplitude ratios are not in a 
best agreement each other for the mass of 14M and mode identification is 
ambiguous between I = and 1 = 1. 

Taken together, we found that the amplitude ratios of ALS 3721 are 
between I = and I = 1 degrees assuming masses in the range of 9M 
and 15M which coincides /3 Cephei instability domain on the HRD. As we 
mentioned above, from the ground based data, we can only determined the 
mode parity of the / degrees in the range of intermediate modes. In general, 
therefore, our results seem in a very good agreement with the theoretical 
expectations. 
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